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Two methods are described for the optimal equalization of a channel 
with "Bump" Equalizers composed of several adjustable-gain Bode Net- 
works. The first method is a general one and applies a steepest descent 
algorithm which minimizes the total mean-squared error (MSE) of the 
equalized channel. It requires continuous gradient information on the 
error-gain relationship in order to determine exactly the optimum equalizer 
adjustments and involves a relatively complicated procedure to calculate 
the gradient. However, the second method, which also applies a steepest 
descent algorithm, develops the necessary gradient information with knowl- 
edge of the error signal only at selected frequencies across the bandwidth 
occupied by the channel. Under idealized assumptions, it is shown that the 
gradients obtained by the second method are exact. When the assumptions 
do not apply exactly, it is shown by computer simulation that the difference 
between the gradients obtained by the two methods is very small. A signifi- 
cant potential advantage of the second method lies in the hardware realiza- 
tion zohich only requires the measurement of the channel error at 2M — 1 
frequencies at the equalizing station (where M is the number of Bode 
Networks in the equalizer). From these frequency domain errors, the 
gradients can be generated as real-time signals and applied to the ap- 
propriate adjustable elements to obtain the optimum gain settings for 
minimum MSE. 

I. INTRODUCTION 

The ideal communication channel exhibits a constant input-output 
gain characteristic over the entire transmission band. In the case of 
a 3/8-inch coaxial cable system, the cable loss varies from 4 dB/mile 
at 1 MHz, to 30 dB/mile at 60 MHz; and to compensate the cable 
loss, repeaters are required at periodic intervals. As is well known, the 
cascaded repeaters cannot exactly compensate the cable loss and this 

1327 



1328 



THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1972 



mismatch results in the so-called "misalignment" of the cable system. 
In addition to the mismatch which is present initially, the channel 
misalignment is affected by the seasonal temperature variation and 
the aging of the components in the system. The objective of the main 
cable equalizers is that, after the equalization, the total input-output 
gain characteristic of the system should be, at all times, as nearly zero 
dB as possible over the entire message band. 

Since the transmission of amplitude information is of major im- 
portance in analog coaxial cable systems, various schemes of input- 
output "amplitude-only" equalization have been studied in the past. 1,2 
It should be noted that the analog signal in the coaxial cable channel 
may contain either voice or digital information, and the transmission 
of voice-type information can be accomplished without phase equaliza- 
tion. For the transmission of high-speed digital information, the neces- 
sary phase equalization is usually furnished in the digital terminals 
and not in the main coaxial cable path. 

In this paper, the Bump Equalizer, 2 which is an "amplitude-only" 
equalizer, is studied and a new adjustment method is presented. The 
Bump Equalizer is composed of a number of Bode Networks, 3 each of 
which can be controlled independently without affecting the other 
networks in the set (Fig. 1). A typical Bode Network is shown in Fig. 2a 
and its transfer function has the form of a bump shape in the frequency 
domain. In the Bump Equalizer of Fig. 3, several Bode Networks are 
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Fig. 1 — Channel equalization with Bump Equalizer. 
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Fig. 2 — Bode Network and its input-output transfer function. 



connected in the feedback and feedforward paths of linear wideband 
amplifiers and the transfer function of the equalizer can be expressed by 

EQL (w) = X 9*B k (w) (dB), (1) 
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Fig. 3 — Equalizer block diagram. 



where w indicates the frequency, g k and B k represent the gain and 
response of the fcth Bode Network respectively, and M is the number 
of Bode Networks in the Bump Equalizer. 

In the past, the so-called "Zero-Forcing" (ZF) method has been used 
for the adjustment of Bump Equalizers. 2 Although the ZF method 
results in zero error at the center frequency of each Bode Network, 
relatively large errors may exist at other frequencies, and as a result 
the ZF method may not be optimal in any overall sense. A better 
error criterion is to minimize the mean-squared error (MSE) over the 
entire bandwidth, and it is this error criterion which is used in this 
paper for equalizer gain adjustment. 

In Section II, the channel and the transfer function of the equalizer 
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are characterized, and several assumptions are made before deriving 
the algorithms which will minimize the MSE of the equalized channel. 
Before the steepest descent method 4 is used to obtain the minimum 
MSE, it is discussed in Section III that the MSE of the equalizer has a 
unique minimum point (no local minima) in the gain parameter space 
of the equalizer, so that this method is assured of finding the true 
minimum MSE at all times. It is shown that under certain conditions, 
the gradient of the MSE, with respect to each gain, can be obtained 
as a real-time signal by measuring the error at only three points in the 
frequency domain. The gradient signal derived can be applied to an 
integrator to adjust the appropriate equalizer gain setting until the 
gradient becomes zero, and the desired optimization is achieved. 
This equalization method is called the "simplified MSE" algorithm. 
Hardware implementation of the algorithm is also discussed in Section 
III. 

Various computer simulations have been carried out and some of 
the results are discussed in Section IV. Both the conventional expres- 
sion for the transfer function of a Bode Network and the measured 
transfer function of a physically realized network have been used in 
the simulation to verify the effectiveness of the derived algorithm 
when used in practical applications. One of the channel misalignments 
used in the simulation resulted from measurements on a working coaxial 
cable system in the field. 

The general steepest descent MSE algorithm is applied to all cases 
to obtain the absolute minimum MSE for each case; and the resulting 
values are compared with the MSE obtained by the simplified MSE 
algorithm. The computer results verify that the simplified algorithm 
derived under the idealized conditions is, in fact, sufficiently close to 
the general algorithm in each case so that the former, which permits 
simplified hardware implementation, can be used as an effective means 
to achieve optimal control of the Bump Equalizer. 

II. CHARACTERIZATION OF CHANNEL AND BUMP EQUALIZER 

The coaxial channel is discussed in this paper principally with respect 
to analog signal transmission. Due to the characteristic of the coaxial 
cable, the bandwidth of the transmitted signal can be quite wide. 
The objective of the equalization discussed is to achieve a constant 
input-output gain characteristic over the entire message band, and 
the transfer function derived is concerned only with the amplitude 
characteristic and not the phase characteristic. Since the transfer 
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function of the Bode Network is symmetric on the "log /" plane, the 
transformed frequency w, to be used in this section is given by log j 
where / is the natural frequency in Hertz. 

2.1 Characterization o] Coaxial Communication Channel 

Let M(iv, t) represent the time-varying channel misalignment which 
is a real-valued function of frequency in units of dB. From the practical 
point of view, however, the channel can be represented as simply M(w), 
since the time-variance can be assumed negligible during the interval of 
any equalization process. Assume also that the Fourier transform of 
the channel misalignment M(w) is limited by a positive constant, 
since the Bump Equalizer to be used is strictly a frequency domain 
equalizer.* Hence, the channel can be characterized in the frequency 
domain by the following series: 

where C„ , j) N , and w„ are certain real numbers, w = log /, and w n+ i — w n 
= l/2pif for all n = 0, 1, • • ■ . 
Eq. (2) also may be expressed as 



M(w) = / X) C« cos (SirMw _ w »)x) dx 

Jo n = 

oo 

^3 C n cos (2irp N w n x) cos (2irp N wx) 

71 = 

4- 2 C- sm (27rpArW„.r) sin (2rp N wx) ? dx 

i-0 > 

\F(x) cos QTrp N wx) + H(x) sin (2irp N wx)} dx, (3) 



'Ml 



l»fl 



where 



F( x ) = 5]C„ cos (2rp N w n x) and H(x) = 2 c '« sin (2TrpA,w n a;). 

n (I »-0 

Since ^ x ^ 1, eq. (3) implies that the shortest frequency domain 
ripple period found in the channel M{w) is l/p N . 



* It should be noted that the Fourier transform of M(w) does not result in an 
impulse response of the channel because of "dB" dimension of M{w). For this 
assumption, however, there is an implicit dual relationship with time domain equal- 
izers, e.g., the tapped delay line equalizers, in which a frequency band limitation of 
the channel is assumed. 
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2.2 Characterization of the Bump Equalizer 

As briefly discussed in the introduction, the frequency domain 
response of the equalizer can be written as 

EQL (w) = £ g k B k (w) (dB), (4) 

where M is the number of Bode Networks in the Bump Equalizer. 
A typical Bode Network is shown in Fig. 2 and its loss is controlled 
by the resistor R. The transfer function, B k (w), can be analytically 
derived and, with a suitable flat-gain amplifier, can be expressed by 
the following equation: 

[E k (l + E k ) + D k (w)Y - D k (w) 
Bk{w) ~ [(l+Etf + DAw)]' (dB) ' (5) 



where 



Rit 



D k (w) = 



(w/w k )H k 



(w/w k ) 2 - 1 ' 
H k = 



C k 



Rc\ k E k 



and 



w k = log(l/2irVllA). 



Since B k {w) of eq. (5) is a quite complicated function of w, the follow- 
ing assumption is made before analyzing the equalizer in detail. 

Assumption 1: Let B k (w) be approximated by 

7T . . IT 



eosmc [- („ - »,)j = — ___ (dB). (6) 

Moreover, if there are M Bode Networks in the equalizer, and if they 
are spaced equally on the iv scale at intervals Aiv* such that Aw = 
w k+l — w k for all k = 1, • • • M — 1, M, then the transfer function of 

* Usually the number of Bode Networks, M, is determined from the practical 
consideration of equalization objectives. p, v defined in eq. (2) determines the Atu 
which is the interval between two adjacent Bode Networks. 
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the equalizer can be expressed by 

EQL (w) = £ g k cosine (^ (to - w k ) J (dB). (7) 

Equations (5) and (6) are plotted in Fig. 2b; and it can be seen that 
cosine (tt/Aiv(w — w k )) approximates the actual transfer function of 
Bode Network given in eq. (5) reasonably well. The maximum dif- 
ference between the two curves of Fig. 2b is 0.0327 dB when \w — w k \ ^ 
Aw and 0.0404 dB when \w — w k \ > Aw. 

III. GAIN OPTIMIZATION USING MEAN-SQUARED ERROR CRITERION 

After a Bump Equalizer has been physically realized and connected 
to the channel, the Aw, which is characteristic of the particular set of 
Bode Networks, cannot be easily altered in the equalizer even though 
the channel misalignment M(w) (and hence p N ) may vary. The optimiza- 
tion here consists of determining the gain parameters g k which will 
minimize the value of MSE defined in this section. One approach to 
the optimization is the employment of the steepest descent method. 4 
In seeking the minimum MSE by this method, the present values of 
g k 's are changed by small amounts in the opposite direction of gradients 
which are the partial derivatives of MSE with respect to each gain 
parameter g k . The process is continued until all the gradients of the 
MSE with respect to the gains g k reach zero or a stationary point. 
Hence, it is implicit in the use of the steepest descent method that the 
surface of MSE in the gain parameter space is a bowl shape, and that 
there exists a unique stationary point which is the global minimum. 
The unique existence of such a stationary point is established before 
a general steepest descent algorithm is derived for the Bump Equalizer; 
then a simplified algorithm is obtained which is shown to be equivalent 
to the general algorithm. Finally, hardware implementation of the 
simplified algorithm is discussed. 

3.1 General Mean- Squared Error Algorithm 

On the dB scale, the residual error after equalization will be 

E(w) = £ g k B k (w) - M(w) (dB). (8) 

If C(w) is the channel characteristic, the channel misalignment is 
defined by M(w) = —C{w) and the MSE can be represented in the 
frequency domain by 
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MSE 



= f E(wf dw. (9) 

V —QQ 



It should be noted that the definition in (9) is a bounded functional, 
since E(w) can be made small as w > w M , where w M is the upper limit 
of the band to be equalized and, moreover, there is complete freedom 
for the assumption of M(w) for w > w M . 

Let G k be a gradient of MSE with respect to gain g k . Then, G k is 
obtained by differentiating MSE with respect to g k : 

G k = ^M§E) = 2 r Bk (w)E(w) dw. (10) 

ogk J-* 

This may be stated as: the gradient G k with respect to gain g k is found 
by cross-correlating the Bode Network function, B k {w), and the error 
function, E{w). The cross-correlation method to obtain the gradient 
in practice has been used elsewhere and can be found in Refs. 5 and 6. 

Theorem 1 (General MSE Algorithm): Let G kJ be the gradient of MSE 
with respect to the kih gain, g k , measured at time t = j (also, let g kii 
indicate the value of the kth gain at time t = j). For the Bump Equalizer, 
the next gain setting of g k (denoted by g k ,j+i) which will reduce the MSE 
is a function of the gain setting g kil and the gradient G kii , as given by 

Qk.i+\ = 9k. i - teQk.i (11) 

for allk = 1,2 ■ • • M where Ac is a small positive constant. As the iterative 
process described by eg. (11) is continued, the gradients G k — > and the 
equalizer reaches the optimum slate. 

Proof: Since an equalizer described by (4) is composed of linearly 
independent networks, there exists a unique set of g k 's which satisfies 
G k = for all k = 1,2, ■ ■ • M, and this set of g k 's results in the mini- 
mum MSE defined in eq. (9) (see Chapter 2 of Ref. 7). Hence, a steepest 
descent algorithm described by (11) must bring the gains to this opti- 
mum stationary point. A general theory on the steepest descent algo- 
rithm is given in Ref. 4. 

Since no specific assumptions were made on the channel and the 
Bump Equalizer, the gradient required for the optimum gain adjust- 
ment may be obtained by eq. (10). 

3.2 Simplified MSE Algorithm and Hardivare Implementation 

In the general steepest descent method described by Theorem 1, 
the gradient G k is obtained by the cross-correlation of the error E(w) 
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and the Bode Network function B k {w). For the hardware implementa- 
tion of the gradient calculation, the error function is multiplied by the 
Bode Network response B k {w) and the product is integrated in the 
frequency domain. Computation of frequency domain cross-correlation 
should be done on a real-time basis, and this fact may prohibit the 
practical application of Theorem 1. However, the following theorem 
shows that the gradient G k can be obtained by measuring E(w) only 
at three different frequencies for each of the given Bode Networks 
B k (w). 

Theorem 2 (Simplified MSE Algorithm): Let the Bump Equalizer satisfy 
Assumption 1 and let the interval Aw between adjacent Bode Networks 
be no greater than the shortest frequency domain ripple period in the 
channel shown in (3), i.e., 

Aw ^ — • (12) 

Pn 

Then, the optimum gain setting of the fcth Bode Network is obtained by 
repeating the following process: 

»,.,♦, = g„, - ac{| *?.(«>. - t) + E < M + \ *«(* + ¥)} • (13) 

where 

k = 1,2,8, ••• M, 

Ac is a small positive constant, and 

Ej(iv k — Aw/ 2), Ej(w h ), and E t (w k + Aw/ 2) are the frequency domain 

errors measured at time t = j at 10 = w k — Aw/2, iv = w h , and w k + 

Aw/2 respectively. 

Proof: The proof is given in the Appendix. 

In the derivation of the results stated in Theorem 2, it is specified 
that the channel has a shortest frequency domain ripple period \/p N 
which is greater than, or equal to, the interval Aw between two adjacent 
Bode Networks. If the channel has ripples having components of 
shorter periods than Aiv, then the gradient shown in Theorem 2 only 
approximates the true gradient which would be obtained by the cross- 
correlation technique in eq. (10). However, the accuracy of the ap- 
proximation depends on the amplitude of ripples whose periods are 
shorter than Aw, and if there exist such ripples of large amplitude, 
one cannot expect a satisfactory equalization even if a general algo- 
rithm is employed. 
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The feedback equation for the optimum gain adjustment using the 
simplified MSE technique is as follows: 

g k (T) = (7,(0) - Ac 

7/6 K* - f ) + em + \ 4» + f )} dt > w 

where k = 1, 2, 3, ■ ■ ■ M, and g k (0) is an initial value of gain g k . The 
hardware implementation of (14) can be achieved in various ways. 
In the block diagram shown in Fig. 3, the input signal source at the 
transmitting station is composed of 2M — 1 unit-amplitude sinusoidal 
oscillators whose frequencies are w t , w t + Aio/2, w 2 , w 2 + Aw/2, • • • , 
w M — Aw/2, and w u . In the block diagram, the band to be equalized 
extends from w = w, to w = w u ■ [It is assumed that E(w x — Aw/2) 
and E(w„ + Aw/2) are zero.] It is possible that in some cases the 
frequencies of the oscillators could be located in the guard bands of 
the channel such that interference between the oscillators and the 
message is avoided and the process could be carried out "in service," 
i.e., in the presence of a "live" message load. At the equalizing station, 
the gradient can be generated by adding the errors with the weighting 
indicated in Theorem 2. Now this gradient, which is a real-time signal, 
is fed to the integrator the output of which can be used to adjust the 
corresponding gain until the optimum condition is reached with respect 
to MSE. In Fig. 3, the function of addition and integration is combined 
by using operational amplifiers. 

TV. RESULTS OF COMPUTER SIMULATION 

In the previous section, two algorithms were derived to obtain the 
MSE optimization. The general algorithm (Theorem 1) can be applied 
for the adjustment of a Bump Equalizer to obtain the minimum MSE, 
but a complex hardware implementation of this scheme may prohibit 
its practical application. Consequently, a simplified algorithm which is 
relatively simple to implement and equivalent to the general one was 
derived. To demonstrate the equivalence of the two algorithms, the 
following two conditions were assumed: 

(i) The channel is represented by a sin x/x series in the frequency 

domain. 
(it) The Bode Network transfer function B k (w) can be approximated 
as in Assumption 1. 

In practice, these restrictions on the channel and the equalizer are 
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Fig. 4 — Channel misalignment, case 1. 



not met exactly. Bode Networks spaced equally in w, for example, 
do not necessarily give the best results. In fact, it has been shown 
experimentally that, when the number of Bode Networks is limited, 
unequally spaced Bode Networks are likely to achieve better equaliza- 
tion. 2 Hence, to verify the effectiveness of the new algorithm, the 
practical limitations were simulated on a digital computer, and the 
resulting performance compared for the two algorithms. 

Case 1: The assumed channel misalignment M(w) is shown in Fig. 4, 
over the natural frequency range from / = 1 MHz to / = 65 MHz. 
Transforming the lower and upper ends of the natural frequency band 
to the logarithmic scale, such that the message band extends from 
w = o (= log e 1) to w = 4.1744 (= log 65), ten Bode Networks are 
specified and spaced equally on the w-axis. The transfer function of 
a physically realized Bode Network was measured and stored in the 
computer for this simulation. 

The results of the simulation are shown in Fig. 5. The total MSE 
of M{w) within the message band is 9.115. The values of MSE obtained 
by the simplified algorithm and the general algorithm are 0.45 and 
0.42, respectively, with the difference resulting mainly from the 
following; 

(t) Since there are sharp corners in M(iv), there exist some ripples 
whose frequency domain periods are shorter than 1/p* [con- 
tradiction to the inequality (12)]. 
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Fig. 5 — Equalized channel, case 1. 

(tt) The measured transfer function of the actual Bode Networks 
used in the simulation differs slightly from the cosine function 
of Assumption 1. 

Case 2: The channel misalignment Af (to) used is one actually measured 
on an existing 20-MHz coaxial cable system, with the bandwidth 
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Fig. 6 — Channel misalignment M(w), case 2. 
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arbitrarily extended to 65 MHz. The misalignment M(iu) used in this 
case is shown in Fig. 6. Ten Bode Networks represented by eq. (5) are 
used in the equalizer, and, in this case, the center frequency of each 
Bode Network is initially optimized for the specified M(w). Con- 
sequently, the resultant array is not spaced equally on the w-axis. 
To apply the simplified MSE algorithm, however, the errors are mea- 
sured at 19 frequencies, 10 of which are the center frequencies of Bode 
Networks and 9 fall between the center frequencies. Total MSE of 
M(w) is 30. Applying the simplified MSE algorithm, MSE = 0.201 
is obtained for the equalized channel and this result is shown in Fig. 7. 
When the general MSE algorithm is used, the absolute minimum 
MSE = 0.186 is obtained. The equalized channel with the general 
MSE algorithm applied is also shown in Fig. 7. 
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Fig. 7 — Equalized channel, case 2. 



V. CONCLUSIONS 

Two algorithms based on the steepest descent method are presented 
in this paper for the optimal gain control of a Bump Equalizer. In 
both cases, the performance index used to evaluate the equalized 
channel is the MSE. The first algorithm is a general MSE algorithm 
and requires MSE gradient information with respect to each gain. 
The required gradients are obtained by a frequency domain cross- 
correlation between the error and the Bode Network to be adjusted. 
For this algorithm to be used, the error signal must be known at all 
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frequencies and this requirement can be a difficult one to implement 
physically. The simplified MSE algorithm derived in Section III, 
however, needs error information at only 2M — 1 frequencies to form 
all of the gradients. Hence, the hardware implementation of the algo- 
rithm is more easily achieved. To derive the algorithm for the simplified 
case, two basic assumptions were made regarding the loss shape of 
the Bode Network and the characteristic of the channel. Under these 
assumptions, the true gradient of the fcth Bode Network is given by 
the weighted sum of the error signal measured at frequencies w k , 
id k + Aw/ 2 , and w k — Aiv/2, where Aid = w k+1 — w k and iv k+l is the 
center frequency of the next higher frequency Bode Network. For the 
hardware realization, the gradient information is applied to integrators, 
the outputs of which in turn control the gain settings, the process being 
continued until all the inputs to the integrator, i.e., gradients, become 
zero (see Fig. 3). 

The computer results given in Section IV show that the reduction 
of total MSE of the equalized channel is negligible by changing from 
the simplified MSE algorithm to the general MSE algorithm, verifying 
the reasonableness of the assumptions made in Section II. 
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APPENDIX 

Proof of Theorem 2 

If we can show that the gradient 

ft = toJ± «(». - f ) + E M + \ E { w „ + f )} , (15) 

the theorem is proved by the result of Theorem 1. From eq. (10), 

G k = 2 J B k (w)E(w) dw 

= 2 f B k (iv)[EQL (w) - M(w)] dw 

J -00 

= 2(G kl - G k2 ), (16) 



where 



G ki = f B k {w) EQL (w) dw 
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and 

G k2 = [ B k (w)M(w) dw. 
By ParservaPs relationship, 8 

G kl = f 2ir&*(0 eql (-t) dt, (17) 

where 

1 r°° 

b k (t) - r- / 5*(w) exp fturt) cfry 

-£(i + «-(f«))-P»M> * M*£ 

= for |«|>i 

and 

1 f" 

eql (0 = o~ / E Q L H ex P O'wO du> 



Ait) 

47T 



(l + cos (^ «jj X) 0< exp O'WiO for | t\ ^ — 



= for U > - — 

1 ' Aw 



Hence, 

2 /.2ir/Au> Jlf 






r 2.r/Au, Jlf / /a ?/ . \\2 

J E flfi exp (j(w k - w<)t)\l + cos [— tjj dt. (18) 



Since w k — w, , = (k — i)6w and the integration with the imaginary 
term in (18) is zero, 



G kl = -£- / £ g, cos ((ft - Aw0 

07T J-2t/Au> t = l 

+ 2 cos (~ M + | cos (A«rf)J d< 

|;2 -2,/a- r 3 i i \ 
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AIOJS , 1 ,1 \ , in , 

= ~2~ \2 9k + 4 g *" 1 + 4 ff * +1 J " ( ^ 

However, 

EQL (u> fc - -yj = - (&_, + ft), EQL (w*) = gr* , 
and 

EQL [w k + -yj = ^ (?* + i + Qk). 
Hence, eq. (19) becomes 

<? M -f [iEQL( Wl -f) + EQL W + i E QL(», + f)]. (20) 
Now 

G 42 = ( 5 t (iy)M(«;) dw. (21) 

From eqs. (3) and (6), 

Gk2 = j cosine f -£- (w — io k )j / {F(x) cos (2rp N wx) 

+ #(.?) sin (27r7? A ri<'.'c) | cte div 

= I cosine f-r— (w — w 4 )J / {/(x) cos (2irp N (w — w k )x) 

+ h{x) sin (2irp s (w — w k )x)\ dx dw, (22) 

where 

F(z) = f(x) cos (2irp N w k ) — h(x) sin (27rp*w*) 

and 

#(z) = j(x) cos (2irp N w k ) + AfcO sin (2irp N iv k ). 

Since cosine (t/Aw(w — if*)) and sin (2irp N {w — w k )x) are even and 
odd functions, respectively, with respect to w = w k axis, 

/ / h(x) cosine (— (w — w k )J sin (2irp N (w — w k )x) div dx = 0. 

Replacing w = u ■+■ u>* and changing "cos" to exponential form, eq. (22) 
becomes 
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- J cosine (~- u) I f(x)(exp (j2irp N ux) + exp (-j2itp N ux)) dx du 

■(exp (j2irp N ux) + exp (— j2irp N ux)) du dx. (23) 
Since 

/ cosine ( — u) exp (jut) du 



and 



Aw ( ' . /Aw .\\ . ■ . | ^ 27r 

= -j- 1 1 + cos l T 1 1 1 for |l|4 = 



= » f » l«l>£, 



eq. (23) becomes 

(?* 2 = -y J /(a)U + cos (ttp^ Awx)} dx for p„.r ^ — 
and 

= for p N x > — ■ 

However, ^ x ^ 1 and Aw ^ l/p^ by definition (12). Hence, 
G k2 = -y / f(*){l + cos (ttp^ Awz)} dx 

- f g Af (». - f ) + WW + 1 lf(<* + f )} ■ (24) 

Combining (21) and (24), the gradient in (16) becomes 
G k = 2(G kl — Gki) 

= A W {|(eQl(»„ - ^) - ^(w, - ^)) + (EQLW - IW) 

+SK(*+f) -«(*+*)} 

- A»g l(« - f ) + E« + | j(« + f )} , (25) 

which is equal to (15). This proves the theorem. 
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